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ABSTRACT 


This  report  presents  the  results  of  a  study  of  the  erosion  of  mate¬ 
rials  by  air-borne  dust.  Erosion  rates  for  several  metals  were  obtained 
over  a  wide  range  of  dust  particle  velocity,  particle  angle  of  impact,  dust 
particle  .size,  and  dust  concentration. 

The  experimental  results  show  that  for  the  range  of  variables  studied, 
erosion  occurs  due  to  the  plastic  displacement  of  material  by  the  dust  par¬ 
ticle.  Relationships  defining  the  effect  of  particle  velocity  and  angle  of 
impact  upon  erosion  rates  are  given.  Experimental  data  show  lug  the  effects  of 
dust  concentration,  dust  size  distribution,  and  material  properties  on  the 
rates  of  erosion  are  presented,  and  the  mechanisms  discussed. 
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L  INTRODUCTION 


1 

Recent  investigations^  at  Southwest  Research  Institute  have  con¬ 
clusively  shown  that  gas  turbine  performance  and  endurance  is  severely 
limited  by  ingestion  of  small  amounts  of  air-borne  dust.  The  results  of 
dust  ingestion  are  erosion  of  shrouds,  turbine  and  compressor  blading,  and 
diffusers  to  the  extent  that  the  machine  may  be  destroyed  in  a  few  hours  of 
operation.  The  trend  towards  use  of  gas  turbine  power-plants  in  military 
ground  vehicles  urgently  requires  a  solution  to  this  problem. 

In  order  to  effect  a  systematic  approach  to  a  solution  to  the  problem 
of  erosion  by  air-borne  dust,  it  became  evident  that  a  basic  understanding 
of  the  mechanisms  of  dust  erosion  was  needed.  To  this  end  this  investiga¬ 
tion  was  directed.  Specifically,  this  investigation  deals  with  the  influence 
of  various  parameters  upon  the  rate  of  erosion  of  several  metals.  These 
parameters  are  dust  particle  velocity,  particle  angle  of  impact,  particle 
size,  particle  concentration  in  the  transporting  air  stream,  and  the  pro¬ 
perties  of  the  material  being  eroded. 

Although  the  primary  objective  of  this  study  is  to  provide  informa¬ 
tion  for  the  improvement  of  gas  turbines  operating  in  a  dusty  atmosphere, 
the  results  can  be  applied  to  a  better  understanding  of  dust  erosion  in 
general,  within  the  range  of  the  various  parameters  used  in  this  investiga¬ 
tion.  Each  parameter  was  varied  in  a  range  as  dictated  by  published  data 
on  the  characteristics  of  natural  air-borne  dust^,  and  on  the  operating 
parameters  of  typical  gas  turbine  machines.  The  ranges  of  the  various 
parameters  in  the  experimental  investigation  were  as  follows; 

_ Parameter _  _ Range _ 

Dust  size  distribution  0-5  micron  to  0-150  micron 

Dust  particle  velocity  150  to  900  fps 

Dust  concentration  0.00027  to  0.0275  gr  dust  /  scf  air 

Dust  particle  angle  of  impact  20  to  90° 

The  eroded  materials  varied  from  lead  (Knoop  hardness  4.6)  to 
martensitic  steel  (Knoop  hardness  750). 


’■'Numbered  superscripts  indicate  correspondingly  numbered  references  in 
List  of  References. 
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II.  PROCEDURE 


A,  Equipment 

The  experimental  setup  is  shown  schematically  in  Figure  1.  The 
air  flow  from  a  220  scfm  air  compressor  is  directed  through  a  shutoff 
valve,  through  a  large  air  dryer  containing  a  desiccant,  through  an  orifice 
plate,  and  then  to  the  air  nozzle.  Upstream  of  the  air  nozzle,  the  required 
weight  rate  of  dust  is  introduced  into  the  air  stream. 

This  flow  of  dust  originates  from  a  dust  feeder,  consisting  of  a 
hydraulic  cylinder,  a  grooved  plate,  and  a  dust  pickup.  Air  pressure  is 
maintained  on  the  piston  end  of  the  cylinder,  and  the  blind  end  of  the  cylin¬ 
der  is  filled  with  water.  Water  is  allowed  to  flow  from  the  blind  end  port 
through  a  needle  valve  and  a  flowmeter,  thence  to  a  drain.  The  speed  of 
piston  travel  may  be  regulated  and  measured  by  the  needle  valve  and  flow¬ 
meter  respectively.  Connected  to  the  end  of  the  piston  is  a  flexible  tube, 
which  leads  to  the  auction  port  of  an  air  ejector.  As  the  piston  retracts,  it 
draws  the  end  of  the  flexible  tube  along  a  dust-filled  groove  in  a  plate  located 
under  the  piston.  The  air  flow  developed  by  the  air  ejector  in  the  flexible 
tube  draws  in  the  dust  from  the  groove,  and  introduces  it  into  the  main  air 
stream.  The  air  ejector  is  powered  by  a  high  pressure  flow  of  air.  The 
pressure  and  flow  rate  of  this  air  flow  are  measured  and  controlled  by  a 
flowmeter,  pressure  gage,  and  a  shutoff  valve. 

The  dust/air  mixture  leaving  the  nozzle  impinges  on  the  material 
specimen  held  in  an  adjustable  specimen  holder.  The  specimen  holder  can 
be  rotated  to  vary  the  angle  of  impact  of  the  dust  particles  on  the  specimen. 
The  axis  of  rotation  lies  in  the  plane  containing  the  specimen  surface  so 
that  the  nozzle  exit- specimen  distance  does  not  vary  with  angle  of  attack. 

The  air -dust  mixture  after  leaving  the  nozzle  is  contained  in  a  shroud,  and 
exhausted  to  atmosphere  outside  the  laboratory. 

B .  Experimental  Procedure 

Prior  to  all  experiments,  it  was  necessary  to  calibrate  the  dust 
feeding  device.  This  was  done  by  filling  the  grooves  in  the  grooved  plate 
with  dry  dust,  and  then  weighing  the  dust  held  in  the  grooves.  After  many 
trials,  a  technique  of  filling  the  grooves  was  found  that  resulted  in  repeat- 
able  dust  weights.  By  dividing  the  dust  weight  by  the  groove  lengths,  a 
value  for  dust  weight  per  unit  length  of  groove  was  obtained.  This  proce¬ 
dure  was  repeated  for  each  dust  size  range  and  dust  type  used  in  the  tests, 


To  Water  Drain 
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The  piston  rate  of  travel  was  correlated  with  the  water  flowmeter 
reading  experimentally.  Then,  by  reading  water  flow  rate  (piston  speed), 
the  weight  rate  of  dust  injected  into  the  main  air  stream  could  be  calculated 
for  any  dust  size  or  type, 

A  typical  run  was  conducted  as  follows;  A  weighed  specimen  was 
installed  in  the  specimen  holder,  and  the  holder  oriented  to  the  desired 
angle.  The  flow  of  high  pressure  air  to  the  ejector  was  started,  and 
adjusted  to  provide  an  adequate  suction  flow  through  the  flexible  tube.  The 
main  shutoff  valve  was  opened  and  adjusted  to  provide  the  required  air 
flow  as  indicated  by  the  orifice  plate  manometers.  Then  the  piston  speed 
was  adjusted  by  manipulation  of  the  needle  valve  to  obtain  the  desired  dust 
feed  rate.  The  test  was  terminated  after  a  certain  weight  of  dust  had  been 
used,  usually  70  or  80  grams.  The  specimen  was  then  reweighed  to  deter¬ 
mine  the  weight  loss.  The  dust  concentration  was  calculated  as 


where  the  symbols  are  as  defined  in  the  List  of  Symbols,  The  total  air 
flow  Q  is  given  as 

Q  =  Ql  +  Q2  (2) 

The  main  air  flow  Qj  was  calculated  using  the  standard  orifice  equations. 
The  dust  injection  air  flow  Q2  was  read  directly  from  the  air  ejector  flow¬ 
meter  and  corrected  to  standard  conditions. 

The  average  air  velocity  at  the  nozzle  exit  was  calculated  directly 
from  the  total  flow  Q,  the  measured  nozzle  area  at  the  exit,  and  the  experi¬ 
mentally  verified  fact  that  the  air  flow  was  essentially  isothermal  downstream 
from  the  main  shutoff  valve.  The  suction  flow  to  the  air  ejector  was  neg¬ 
lected,  since  it  was  found  to  be  small  compared  to  the  measured  flows. 

A  survey  of  the  air  velocity  at  the  nozzle  exit  was  made  with  a 
pitot  tube.  This  served  to  verify  the  calculated  average  air  velocity,  and 
showed  that  the  air  velocity  was  essentially  constant  over  the  extent  of  the 
specimen.  From  the  pitot  tube  data,  a  factor  was  obtained  to  convert 
average  air  velocity  experienced  by  the  specimen;  in  other  words,  a  cor¬ 
rection  for  nozzle  boundary  layer  effects  was  used. 

Specimen  weight  losses  due  to  erosion  are  presented  as  specimen 
weight  loss  per  unit  specimen  area,  per  unit  weight  of  dust  impinging  on 
the  specimen.  This  factor  is  denoted  by  the  term  "weight  loss  factor.  " 

The  "volume  loss  factor"  is  the  weight  loss  factor  divided  by  the  specimen 
material  density. 
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C.  Test  Specimens 

Table  I  presents  the  properties  o£  the  materials  tested. 

Prior  to  each  test,  the  specimen  was  polished  to  remove  surface 
scratches,  and  then  weighed  on  an  analytical  balance.  The  specimen  was 
weighed  twice  before  the  test,  and  twice  afterwards,  the  weight  loss  being 
computed  from  the  average  weights.  Weights  were  recorded  to  the  nearest 
tenth  of  a  milligram,  and  the  average  weight  seldom  differed  from  the 
extreme  by  more  than  two  tenths  of  a  milligram. 

The  specimens  were  stored  in  a  desiccator,  and  were  handled  by 
tweezers  to  prevent  contamination  of  the  surfaces.  Before  weighing,  the 
specimens  were  carefully  cleaned  to  remove  air -borne  particles. 

Each  specimen  was  approximately  3/8  in.  square,  and  .060  in. 
thick.  This  relatively  small  size  was  dictated  by  considerations  of  nozzle 
area,  air  velocity,  and  available  compressor  capacity.  The  resulting 
small  specimen  weight  contributes  to  more  accurate  weight  loss  measure¬ 
ments,  since  for  heavy  specimens  the  weight  loss  is  the  small  difference 
between  two  large  numbers,  with  the  attending  inaccuracies, 

D.  Dust 


For  the  majority  of  the  erosion  experiments,  the  dust  used  was 
silica  flour,  with  the  following  typical  chemical  analysis: 


Component _  %  by  Weight 


Silicon  dioxide  (Sj^02)  99.90 

Iron  Oxide  (Fe203)  0,018 

Aluminum  Oxide  {AI2O3)  0.12 

Titanium  dioxide  (Ti02)  0,007 

Calcium.  Oxide  (CaO)  0.01 

Magnesium  Oxide  (MgO)  Trace 


The  silicon  dioxide  in  this  dust  is  present  in  the  form  of  quartz, 
with  a  Moh's  hardness  of  7. 

For  comparison  purposes,  several  tests  were  run  using  the  two 
grades  of  SAE  test  dust^  with  the  following  chemical  analysis: 


6 


TABLE  I.  SPECIMEN  MATERIAL  PROPERTIES 


Material 

Knoop 
Hardness, 
KG/ mm2 

Density, 

gr/cc 

Tensile 
Strength, 
psi  X  10"2 

Elongation 
in  2  Inches 
% 

C-1050  Steel 
(Martensitic) 

750* 

7.  80 

375 

5 

C-1050  Steel 
(Tempered) 

336* 

7.80 

128 

23 

C-1050  Steel 
(Pearlitic) 

192* 

7.80 

102 

30 

303  Stainless  Steel 

190* 

7.93 

100 

50 

Copper 

63* 

8.92 

34 

45 

6061 -T6  Aluminum 

65* 

2.65 

45 

12 

6061-0  Aluminum 

43* 

2.  65 

18 

25 

Lead 

4.  6* 

11.3 

1.9 

100 

’•■■Measured  values  - 

all  other  values 

estimated  from  handbook  data. 

(All  Knoop  hardness  numbers  determined  using  200  gram  load.) 
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Component 


% 


by  Weight 


Silicon  dioxide  (SiO^)  67  to  69 

Iron  Oxide  (Fe203)  3  to  5 

Aluminum  Oxide  (AI2O3)  15  to  17 

Calcium  Oxide  (CaO)  2  to  4 

Magnesium  Oxide  (MgO)  0.5  to  1.5 

Total  Alkalis  3  to  5 

Ignition  Loss  2  to  3 


More  significant  in  dust  erosion  studies  is  the  mineral  identity  of 
the  dust,  as  shown  below  for  the  SAE  test  dust: 


Mineral  Identity 


%  by  Weight  Moh' s  Hardness 


Plagioclase  Feldspars  65  -  70  6  -  6,5 

Quartz  15-18  7 

Amphibole  Hornblende  10-12  5-6 

Vesuvianite  0.5-1  6.5 

Muscovite  Mica  Approx.  0.5  2,5-3 

Biotite  Mica  Approx.  0,5  2,5-3 

Tourmaline  Approx.  0.5  7-7,5 


The  size  distributions  of  the  dusts  used  in  this  study  are  shown  in 
Figure  2. 


E.  Accuracies  of  Controlled  Parameters 


To  determine  the  effect  of  the  various  parameters  upon  the  rate 
of  erosion,  the  parameter  being  investigated  was  assigned  a  different  value 
for  each  test  in  a  series  of  testa  while  the  other  parameters  were  held  at 
a  constant  value  for  the  complete  series.  The  dependent  variable,  weight 
loss,  was  measured  for  each  test.  Certain  inherent  characteristics  of 
the  experimental  setup  introduced  a  margin  of  error  in  the  value  of  these 
controlled  parameters. 

The  primary  source  of  error  in  the  controlled  parameters  arises 
from  the  dynamics  of  the  dust  particles  in  the  air  stream.  The  particles 
are  introduced  into  the  stream  with  some  initial  velocity  and  are  then 
accelerated  by  the  air  stream  until  they  impinge  upon  the  specimen.  The 
impinging  velocity  of  the  particles  depends  upon  their  size,  weight,  shape, 
and  upon  the  velocity,  pressure,  and  temperature  of  the  airstream.  Other 
factors  being  constant,  the  smaller  particle  will  more  nearly  attain  the 
airstream  velocity  than  the  larger  particle. 
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In  the  neighborhood  of  the  specimen,  the  airstream deflects  around 
the  specimen.  Dust  particles  in  this  stream  are  thus  subjected  to  drag 
forces  acting  to  change  the  particle  direction  of  travel.  This  causes  the 
dust  particles  to  strike  the  specimen  at  angles  other  than  the  test  angle. 
Throughout  this  report,  the  test  impact  angle  will  be  defined  as  the  angle 
between  the  specimen  and  the  nozzle  centerline..  The  actual  angle  of 
impact  is  defined  as  the  angle  between  the  velocity  vector  of  a  particular 
particle  and  the  specimen  at  the  instant  of  contact.  The  difference  between 
these  two  angles  increases  with  decreasing  particle  size,  decreasing  air 
velocities,  and  increasing  impact  angles. 

An  estimation  of  the  velocity  error  may  be  obtained  from  a  solution 
of  the  particle  one-dimensional  equation  of  motion.  From  Dallavalle^,  the 
drag  coefficient  for  spherical  particles  is 

40 

Cr  =  0.  4  +  ^  (2  <  R  <  500)  (3) 

R 


where 


Po*^ 

R  = - (W  -  V)  (4) 

It  may  be  shown  (see  Appendix  A)  that  for  dust  particles  less  than 
approximately  150  microns  in  diameter  in  an  airstream  of  normal  tempera¬ 
ture,  the  distance -velocity  relationship  of  the  dust  particle  is  given  by 


Pq  X 
P  d 


30 


R  4  100  ^ol  ^  ,  1  ,  /  R  +  lOO  \ 
Ro  +  100  R  J  ^  3  \  Ro  +  100  / 


(5) 


where 


R 


(6) 


The  drag  coefficient  for  higher  particle  Reynold's  numbers  is 
given'^  by 

Cr  =  0.  44  (500  <  R  <  10^)  (7) 


Thus,  it  is  seen  that  Equation  (5)  is  more  approximate  but  still  useful  at 
Reynold's  numbers  greater  than  the  range  indicated  by  Equation  (3),  Solu¬ 
tions  to  Equation  (5)  are  shown  graphically  in  Figure  3,  for  the  nozzle 


Particle  Velocity  at  Nozzle  Exit  -  fps 


length  used  in  the  experiments  (5  feet)  and  for  various  particle  sizes.  The 
consequences  of  the  deviation  of  particle  velocity  from  the  nozzle  air 
velocity  will  be  considered  elsewhere  in  this  report. 

The  determination  of  the  true  impact  angle  of  the  dust  particle  on 
the  specimen  presents  considerable  difficulties.  An  experimental  technique 
was  devised  to  make  a  rough  approximation  of  the  maximum  error  in  impact 
angle.  A  glass  specimen  was  placed  in  the  specimen  holder,  oriented  to 
obtain  a  test  impact  angle  of  90“ ,  and  impinged  with  a  small  amount  of  0-5 
micron  dust  at  approximately  500  feet  per  second  air  velocity.  The  glass 
specimen  was  then  inspected  under  a  microscope,  and  the  number  and 
location  of  the  resulting  pits  in  the  glass  recorded.  For  a  first  approxi¬ 
mation,  it  can  be  shown  that  (see  Appendix  B) 

N 

sin  a  =  -  sin  6  (8) 

^max 

Using  Equation  (8),  a  was  calculated  to  be  approximately  40°,  the  error  in 
this  case  being  roughly  60“. 

For  dust  in  the  size  range  50-74  microns,  at  the  same  air  velocity 
and  test  angle,  a  was  found  to  be  essentially  90";  that  is,  the  actual  angle 
was  approximately  equal  to  the  test  angle. 

The  60“  error  for  the  0-5  micron  dust  is  considered  an  extreme 
case,  and  the  average  error  for  the  large  majority  of  the  tests  is  probably 
much  lower  than  this  value.  The  effects  of  this  error  can  be  seen,  however, 
and  will  be  considered  again  in  the  discussion  of  the  test  results. 
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in.  RESULTS 


A .  Th  eoretic  al 


In  a  paper  by  Finnie®,  a  theoretical  analysis  of  the  volume  of  mate¬ 
rial  removered  by  an  impinging  abrasive  grain  is  presented.  A  summary  of 
this  analysis  is  presented  below. 

Figure  4  shows  an  idealized  representation  of  an  abrasive  particle 
removing  material  from  a  surface.  The  objective  of  the  analysis  is  to 
determine  the  equation  of  motion  of  the  particle  through  the  surface,  thus 
allowing  an  estimation  of  the  volume  of  material  removed. 


V 


FIGURE  4.  REMOVAL  OF  MATERIAL  BY 
AN  ABRASIVE  PARTICLE 

Several  assumptions  are  required.  First,  it  is  assumed  that  the 
particle  has  no  initial  rotation.  It  is  very  probable  that  actually  the  particles 
have  rotations  in  all  directions,  but  it  seems  reasonable  to  take  the  assump¬ 
tion  as  an  average  condition. 

Since  it  may  be  shown  that  the  particle  rotation  during  cutting  is 
small,  it  is  assumed  that 


Yt  “  Vq 

xt  *  34q  +  r<b 


(9) 


’3 


It  is  further  asBumed  that  the  force  on  the  particle  during  cutting  will  be 
proportional  to  the  area  of  contact  times  a  plastic  flow  stress  p.  The  ratio 
of  the  vertical  to  hori/iontal  force  on  the  particle  (K)  is  assumed  to  be  con¬ 
stant  during  cutting,  in  view  of  the  first  and  second  assumptions. 

Then,  the  equations  of  motions  may  be  written  as 

my-p  +  pKF(yp')  =  0  (10) 


mx-p  +  p 


(l  t2^)F(yT)  =  0 


(11) 


The  term  F(y,p)  is  the  projected  area  of  particle-surface  contact  on  the 
vertical  plane.  The  initial  conditions  are  given  as 

x-p  =  V  cos  a  (12) 

y-p  =  V  sin  a  (13) 

The  incremental  volume  removed  is 


F(yT) 

dq  =  —  dxrp  (H) 

where  i(j  is  the  ratio  of  the  projected  area  of  contact  to  the  projected  area  of 
contact  underneath  the  original  surface. 

Integrating  Equation  (14)  after  substituting  for  F(y'p)  from  Equation 
(1 1)  yields 


(15) 


where  xip  is  the  horizontal  velocity  of  the  particle  tip  when  cutting  is  complete. 


One  end  condition  is  xp  =  0,  then 


q  = 


mV^ 


cos^  a 


(16) 
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If  the  particle  tip  still  has  horizontal  velocity  when  cutting  is  complete,  it 
can  be  shown  that  for  this  case 


q  = 


2  .  _  2  sin^  a 

p  [sin  2a - - - 


where 


P  = 


K 


1  + 


mr‘ 


(17) 


(18) 


Equations  (16)  and  (17)  are  equal  for  tan  a  =  P/2.  Maximum  volume  removal 
occurs  for  tan  2a  =  P.  Equation  (16)  is  valid  for  tan"^  P/2  <  a  <  90"  and 
Equation  (17)  is  to  be  used  for  0  <  q<  tan“^  P/2.  Equations  (16)  and  (17) 
predict  that  for  a  single  particle,  the  volume  removed  is  proportional  to  the 
particle  velocity  squared,  and  inversely  proportional  to  the  flow  stress  of 
the  material,  In  addition,  it  is  predicted  that  the  volume  removed  will 
exhibit  a  maximum  when  plotted  versus  impact  angle  a,  the  angle  for  maxi¬ 
mum  erosion  dependent  on  the  parameter  P,  Equation  (17)  shows  that  for 
an  impact  angle  of  zero  degrees,  the  volume  loss  is  zero.  Also  for  an 
impact  angle  of  90®,  Equation  (16)  predicts  no  volume  loss.  Figure  5  shows 
the  theoretical  results  plotted  for  several  values  of  the  parameter  P, 
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FIGURE  5,  PREDICTED  VOLUME  REMOVAL,  FOR 
A  SINGLE  ABRASIVE  PARTICLE 
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B.  Angle  of  Impact  Effect 

The  effect  of  angle  of  impact  upon  material  erosion  will  be  presented 
using  the  weight  loss  factor  as  a  measure  of  the  erosion  rate.  Before  dis¬ 
cussing  the  results  obtained,  the  method  of  calculating  this  factor  should  be 
considered.  As  previously  explained,  the  weight  loss  factor  is  the  weight 
of  material  eroded  per  unit  specimen  area  per  unit  weight  of  dust  impinged 
on  the  specimen.  From  Figure  6  it  is  seen  that  the  weight  of  dust  impinging 
on  the  specimen  varies  with  the  angle  of  impact. 


FIGURE  6.  NOZZLE-SPECIMEN  GEOMETRY 

The  weight  of  dust  striking  the  specimen  is  given  by  (total  weight  of  dust 

(A,  sin  0  \ 

- j  .  Then  the  weight  loss  factor  is  expressed  as 

_ weight  eroded _  _  (weight  eroded)  (a) 

(A)  (wt  of  dust  striking  specimen)  (A^)  (total  dust  wt)  (sin  6) 

The  effect  of  test  angle  of  impact  upon  erosion  rate  is  shown  in 
Figures  7  through  15.  In  several  of  the  graphs,  there  is  considerable 
scatter  of  data,  especially  at  low  angles  of  impact.  The  manner  in  which 
to  draw  the  curves  to  best  fit  these  data  is  certainly  not  obvious,  Howei/er 
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Test  Impact  Angle  -  Degrees 


FIGURE  7.  WEIGHT  LOSS  FACTOR  VERSUS  TEST  IMPACT 
ANCLE,  MARTENSITIC  C-1050  STEEL,  AIR 
VELOCITY  475  FPS,  DUST  CONCENTRATION 
0.011  GR/Ft3,  0  -  74  MICRON  SILICA  FLOUR 
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^  ^  Test  Impact  Angle  -  Degrees 

FIGURE  8.  WEIGHT  LOSS  FACTOR  VERSUS  TEST  IMPACT 
ANGLE,  PEARLITIC  C-1050  STEEL,  AIR 
VELOCITY  475  FPS,  DUST  CONCENTRATION 
0.011  GR/FT3,  0-74  MICRON  SILICA  FLOUR 
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Test  Impact  Angle  -  Degrees 


FIGURE  9.  WEIGHT  LOSS  FACTOR  VERSUS  TEST  IMPACT 
ANGLE,  TEMPERED  C-1 050  STEEL,  AIR 
VELOCITY  475  FPS,  DUST  CONCENTRATION 
0,011  GR/Ft3,  0-74  MICRON  SILICA  FLOUR 
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FIGURE  10,  WEIGHT  LOSS  FACTOR  VERSUS  TEST  IMPACT 
ANGLE,  MARTENSITIC  C-1 050  STEEL,  AIR 
VELOCITY  835  FPS,  DUST  CONCENTRATION 
0.013  GR/Ft3,  0  -  74  MICRON  SILICA  FLOUR 
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FIGURE  11.  WEIGHT  LOSS  FACTOR  VERSUS  TEST  IMPACT 
ANGLE,  PEARLITIC  C-1050  STEEL,  AIR 
VELOCITY  835  FPS,  DUST  CONCENTRATION 
0.013  GR/FT^,  0-74  MICRON  SILICA  FLOUR 
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FIGURE  12.  WEIGHT  LOSS  FACTOR  VERSUS  TEST  IMPACT 
ANGLE,  TEMPERED  C-1050  STEEL,  AIR 
VELOCITY  835  FPS,  DUST  CONCENTRATION 
0.013  GR/Ft3,  0-74  MICRON  SILICA  FLOUR 
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FIGURE  13.  WEIGHT  LOSS  FACTOR  VERSUS  TEST  IMPACT 
ANGLE,  6061 -T6  ALUMINUM  ALLOY,  AIR 
VELOCITY  835  FPS,  DUST  CONCENTRATION 
0.013  GR/FT3,  0-74  MICRON  SILICA  FLOUR 
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FIGURE  14.  WEIGHT  LOSS  FACTOR  VERSUS  TEST  IMPACT 
ANGLE,  PEARLITIC  C-1050  STEEL,  AIR 
VELOCITY  835  FPS,  DUST  CONCENTRATION 
0.013  GR/FT^,  SAE  TEST  DUST.  COARSE  GRADE 
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FIGURE  15.  WEIGHT  LOSS  FACTOR  VERSUS  TEST  IMPACT  ANGLE, 
PEARLITIC  C-1050  STEEL,  AIR  VELOCITY  835  FPS, 
DUST  CONCENTRATION  0.013  GR/FT^,  SAE  TEST 
DUST,  FINE  GRADE 
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from  the  expression  for  the  weight  loss  factor,  it  is  seen  that  the  sin  0 
appears  in  the  denominator.  At  low  impact  angles,  a  small  change  in  the 
test  angle  of  impact  6  has  a  large  effect  upon  the  reciprocal  of  the  sin  0. 

Thus,  a  small  error  in  test  angle  of  impact  causes  a  small  error  in  the 
location  of  the  datum  point  along  the  abscissa,  and  a  relatively  large  error 
in  its  location  along  the  ordinate.  Small  errors  in  test  angle  were  difficult 
to  prevent,  the  small  size  of  the  specimen  making  very  accurate  installation 
and  alignment  difficult.  Erosion  of  the  specimen  holder  also  contributed  to 
errors  in  test  angle.  To  aid  in  plotting  the  data,  the  weight  loss  factor  for 
each  point  was  multipled  by  the  sin  of  the  corresponding  impact  angle,  and 
these  points  were  plotted  with  dashed  lines.  A  best  fit  curve  was  easily 
drawn  through  these  points,  since  errors  in  the  weight  loss  factor  due  to 
test  impact  angle  errors  were  greatly  decreased.  The  points  on  the  dashed 
curve  were  then  divided  by  sin  9  and  thus  transferred  to  weight  loss  factor 
coordinates  as  the  curve  of  best  fit. 

Figures  7  through  15  show  the  trends  predicted  by  Finnic' s  deriva¬ 
tion.  As  the  angle  of  impact  decreases  from  90",  the  erosion  rate  increases 
to  a  maximum,  then  abruptly  drops  towards  zero  at  an  impact  angle  of  zero 
degrees.  Agreement  between  the  theoretical  analysis  and  experiment  is 
generally  good,  the  poorest  agreement  occurring  at  high  impact  angles.  This 
is  due  to  the  fact  that  at  high  test  impact  angles,  the  dust  particles  strike 
the  specimen  at  angles  lower  than  the  test  angle  due  to  the  change  in  the 
direction  of  the  air  stream  in  the  neighborhood  of  the  specimen.  In  addition, 
the  surface  of  the  specimen  is  pitted  from  prior  impacts,  and  the  actual 
impact  angle  of  a  given  particle  will  be  affected  by  the  shape  of  these  pits. 
These  two  factors  contribute  to  the  excess  of  experimental  erosion  over  the 
predicted  trend. 

The  data  shown  represent  tests  conducted  using  two  different  air- 
stream  velocities.  The  data  taken  at  the  lower  velocity  (Figs.  7-9) 
exhibit  maximum  erosion  rates  at  an  erosion  angle  of  approximately  20". 

The  higher  air  velocity  plots  (Figs.  9  -  15),  on  the  other  hand,  indicate 
that  maximum  erosion  occurs  at  an  angle  in  the  neighborhood  of  10". 

From  Equation  (18)  and  the  subsequent  discussion,  it  is  seen  that  the  angle 
of  maximum  erosion  is  a  function  of  the  geometry  of  the  eroding  particle, 
and  does  not  depend  on  the  impact  velocity.  It  would  be  expected,  there¬ 
fore,  that  the  angle  for  maximum  erosion  would  be  equal  for  the  two  types 
of  tests.  The  reason  for  this  inconsistency  seems  to  be,  again,  a  conse  ¬ 
quence  of  the  actual  impact  angle  being  consistently  smaller  than  the  test 
angle.  For  lower  velocities,  this  effect  is  more  pronounced,  and  results 
in  the  maximum  point  being  shifted  towards  higher  angles.  Thus,  it  seems 
probable  that  the  correct  angle  for  maximum  erosion  is  more  nearly  that 
shown  for  the  higher  velocity  tests  than  for  the  lower. 
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Figures  14  and  15  show  data  obtained  using  SAE  test  dust,  coarse 
grade  and  fine  grade  respectively.  All  the  other  curves  represent  data 
resulting  from  erosion  by  the  silica  flour.  It  is  seen  that  the  general  shape 
of  the  curves  for  SAE  test  dust  is  somewhat  different  from  those  for  silica 
flour,  although  the  angles  for  maximum  erosion  rate  are  approximately 
equal.  This  is  construed  to  mean  that  the  values  of  P  [see  Eq.  (18)]  for 
each  dust  type  are  approximately  the  same,  since  the  angle  of  maximum 
erosion  depends  only  on  the  value  of  P.  Since  the  SAE  test  dust  is  a  natural 
dust  that  has  been  subjected  to  weathering,  and  the  silica  flour  is  artificially 
crushed,  it  would  be  expected  that  the  silica  flour  would  exhibit  a  higher 
degree  of  angularity  than  the  SAE  test  dust.  Microscopic  examination  does 
not  reveal  a  marked  difference  in  angularity  between  the  two;  however, 
undetected  differences  in  angularity  may  contribute  to  significant  differences 
in  erosion  capabilities.  In  addition,  mineral  composition  and  particle  size 
distribution  differ  in  the  two  dust  types,  and  this  doubtlessly  contributes  to 
the  difference  in  the  shape  of  the  two  sets  of  curves. 

C.  Erosion  Wear  Patterns 


An  interesting  phenomenon  observed  during  dust  erosion  tests  on  gas 
turbines  previously  conducted  at  Southwest  Research  Institute^  was  the  sur¬ 
face  ripples  formed  on  gas  turbine  components  that  had  experienced  high 
rates  of  erosion.  These  ripples  were  particularly  striking  on  the  compressor 
impeller  (2025T6  aluminum  alloy)  and  compressor  diffuser  (QQ-A-601), 

Class  3M,  Condition  T6  Aluminum  Alloy  .  The  ripples  were  observed  to 
form  in  locations  where  the  impact  angle  was  low  and  the  particle  velocity 
high.  The  ripples  extend  in  roughly  parallel  lines  perpendicular  to  the 
direction  of  air  flow.  No  ripples  were  observed  with  small  size  ranges  of 
dust,  but  were  common  when  the  dust  used  was  in  the  size  range  of  0  -  74 
microns. 

Similar  ripple  patterns  were  observed  during  the  present  series  of 
tests.  The  first  observation  was  made  on  martensitic  C-1050  steel  speci¬ 
mens  subjected  to  erosion  at  a  90°  impact  angle.  The  patterns  were  rough 
concentric  circles  with  the  center  of  the  circles  coinciding  with  the  speci  ¬ 
men  center.  The  amplitude  of  the  ripples  were  on  the  order  of  10  microns, 
the  wavelength  approximately  4  mm.  After  conducting  erosion  tests  on  a 
lead  specimen  at  an  impact  angle  of  40°,  relatively  deep,  well-defined 
patterns  were  visible.  Figure  16  is  a  photograph  of  this  specimen.  The 
ripples  represented  in  Figure  16  have  a  wavelength  of  approximately  0.  16  mm. 

Consideration  of  the  equations  of  motion  for  a  single  abrasive  particle 
leads  to  a  possible  explanation  of  the  formulation  of  these  wear  patterns. 

From  the  equations  of  motion  for  a  single  abrasive  particle  (Eqs.  9  and  10). 


FIGURE  16.  EROSION  RIPPLE  PATTERN-LEAD  SPECIMEN 
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the  horizontal  and  vertical  components  of  the  particle  velocity  may  be 
found  at  the  instant  the  particle  leaves  the  specimen  surface 

I  Z 

i-p  =  V  cos  a  -  —  V  sin  a 

y,p  =  -V  sin  a 

These  results  are  valid  for  0  <  a<  tan"  ^  P/2,  For  small  values  of  a,  it 
is  seen  that  the  horizontal  velocity  component  may  take  on  significant 
values.  It  is  thought  that  the  ripple  wear  patterns  are  a  result  of  the  skipping 
motion  of  the  individual  particles,  a  motion  much  like  that  of  a  flat  stone 
skipping  across  a  calm  surface  of  water.  For  the  case  of  the  dust  particle, 
however,  the  particle  will  gain  additional  velocity  from  the  air  stream  as 
it  rises  after  leaving  the  specimen  surface.  The  particle  is  then  directed 
back  down  onto  the  specimen  surface  by  the  airstream.  In  experiments 
with  lead  specimens  oriented  at  40“  to  the  airstream  the  ripples  were 
found  to  be  deepest  and  most  distinct  adjacent  to  and  downstream  from  a 
discontinuity  consisting  of  a  gap  between  the  specimen  and  the  specimen 
holder.  Deep  scratches  cut  into  into  the  lead  specimen  caused  similar 
deep  distinct  ripples  leading  to  the  theory  that  the  discontinuities  are  the 
origin  of  the  ripple  patterns.  Once  the  pattern  is  originated,  it  is  stable, 
since  each  ripple  becomes  a  source  for  further  rippling  downstream.  The 
rippling  pattern  is  observed  to  decay  and  lose  its  regularity  at  points  far 
from  the  originating  discontinuity.  As  erosion  continues,  the  ripples 
become  deeper  and  more  distinct. 

At  first,  the  concentric  ripples  on  the  martensitic  steel  specimen 
did  not  seem  to  reinforce  this  theory  due  to  the  absence  of  an  originating 
discontinuity.,  However,  microscopic  examination  of  several  martensitic 
steel  specimens  showed  that  all  had  a  very  shallow  dimple  in  the  center  of 
the  specimen,  probably  a  result  of  the  severe  heat  treating  quench.  Thus, 
particles  striking  near  the  center  of  the  specimen  were  then  accelerated 
along  the  surface  of  the  specimen  by  the  airstream,  encountered  the  edge 
of  the  dimple,  and  rebounded  to  begin  the  wear  pattern. 

It  should  be  noted  that  the  presence  of  a  deep  ripple  pattern  on  the 
specimen  surface  would  be  expected  to  confuse  the  effect  of  test  impact 
angle  upon  erosion  rates.  This  was  not  verified  since  no  angle  of  impact 
data  was  obtained  for  lead,  and  the  depth  of  the  ripple  pattern  on  other 
specimens  was  so  small  as  to  have  negligible  effect  on  erosion  rates. 

The  fact  that  each  dust  particle  makes  multiple  impacts  on  the 
specimen  surface  will  tend  to  make  the  test  angle  of  impact  versus  erosion 
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rate  curves  more  flat,  and  will  contribute  to  the  observed  erosion  at  a  test 
angle  of  90°  ,  in  contradiction  to  the  theoretical  analysis. 

D,  Velocity  Effect 

From  the  derivation  for  material  removal  by  a  single  abrasive 
particle,  it  would  be  expected  that  the  material  removed  will  vary  as  the 
particle  velocity  squared.  The  results  of  several  erosion  tests  with  air 
velocity  as  the  independent  variable  are  shown  in  Figure  17  through  20. 
These  results  show  that  the  predicted  relationship  is  borne  out  extremely 
well  for  a  test  impact  angle  of  40°,  and  less  well  for  a  test  impact  angle  of 
90°. 


At  first  consideration,  it  would  seem  that  the  differences  in  air 
velocity  and  actual  particle  velocity  would  cause  the  results  shown  to  be 
in  considerable  error  when  considering  the  effects  of  actual  particle 
velocity  upon  erosion  rate.  However,  Figure  3  shows  that,  to  a  very  good 
approximation,  the  particle  velocity  is  a  linear  function  of  the  airstream 
velocity.  Denoting  the  various  sizes  of  particles  in  the  airstream  by 
subscripts 


Vi  =  C]  W 

V2  =  C2W 


Va  =  C^W 


(]9) 


where  C  is  a  constant  of  proportionality.  Assuming  that  for  each  size  of 
particle,  the  erosion  rate  varies  as  the  nth  power  of  the  velocity,  the 
weight  loss  factor  is  given  by 


.  + 


(20) 


where  k  is  a  constant.  Substituting  Equation  (19)  into  Equation  (20)  and 
collecting  terms  gives 


w  =  (kiC^mi  +  k^C^m:: 


n 


(21) 


or  w  =  (constant)  (W^),  This  shows  that  the  slope  of  the  weight  loss -air- 
stream  velocity  curve  will  accurately  reflect  the  actual  weight  loss- 
particle  velocity  power  relationship.  The  difference  in  airstream  velocity 
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FIGURE  18.  WEIGHT  LOSS  FACTOR  VERSUS  AIR  VELOCITY 

PEARLITTC  C-1050  STEEL,  DUST  CONCENTRATION 
0.  013  GR/FT^,  0  74  MICRON  SILICA  FLOUR 
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FIGURE  20.  WEIGHT  LOSS  FACTOR  VERSUS  AIR  VELOCITY 

6061-T6  ALUMINUM  ALLOY,  DUST  CONCENTRATION 
0.  013  GR/Ft3,  0  -  74  MICRON  SILICA  FLOUR 


30 


and  particle  velocity  will  not,  therefore,  cause  an  appreciable  error  in  the 
slope  of  the  weight  loss-air  velocity  curves. 

The  difference  in  slope  of  the  weight  loss-air  velocity  curves  for  a 
given  material  at  test  angles  of  40°  and  90°  is  believed  to  be  primarily  due 
to  the  difference  in  test  angle  and  particle  impact  angle.  This  error,  pre¬ 
viously  discussed,  causes  the  data  for  the  90°  test  angle  to  actually  reflect 
the  results  of  impacts  varying  from  90°  for  the  larger  particles  to  some 
lower  angle  for  the  smaller  particle.  This  difference  in  test  angle  and 
actual  impact  angle  changes  with  air  velocity,  and  since  the  weight  loss 
varies  with  impact  angle  for  a  given  particle,  the  results  shown  for  a  90° 
test  angle  probably  do  not  accurately  represent  the  true  velocity  effect. 

This  error  is  greatly  reduced  for  lower  test  angles,  and  the  data  for  the 
40°  test  angle  provide  a  more  accurate  representation  of  the  relationship. 

The  above  results,  with  the  results  of  the  effect  of  variation  in 
angle  of  impact,  in  general  show  good  agreement  with  the  equations  derived 
by  Finnie.  These  equations  describe  an  erosion  process  characterized  by 
plastic  deformation  of  the  material.  In  another  paper^,  Finnie  has  described 
the  characteristics  of  erosion  of  brittle  materials.  The  brittle  materials 
exhibit  maximum  erosion  at  impact  angles  of  approximately  90°,  the  erosion 
rate  then  decreasing  monotonously  to  zero  at  a  0°  impact  angle.  In  addition, 
the  erosion  rate  of  brittle  materials  varies  as  the  particle  velocity  taken  to 
a  power  greater  than  two.  For  instance,  it  was  found  that  erosion  of  porcelain 
showed  w  ■■  and  in  tests  of  brittle  steel,  the  exponent  varied  between 

5.  5  and  6.  5.  In  the  present  study,  the  brittle  type  of  erosion  was  not 
encountered.  It  is  interesting  to  note  that  during  erosion  tests  conducted  by 
Finnie  using  particles  of  60-mesh  (250  micron)  silicon  carbide,  it  was 
found  that  C-1055  SAE  steel  hardened  to  Rockwell  C  63  showed  the  charac¬ 
teristics  of  a  brittle  material.  In  the  present  study,  a  nearly  identical 
material  (martensitic  C-1050  SAE  steel,  Rockwell  C  62)  shows  the  erosion 
characteristics  of  a  ductile  material.  This  indicates  that  the  line  of  separa¬ 
tion  between  brittle  and  ductile  materials  is  not  a  definite  one,  and  may 
depend  a  great  deal  on  abrasive  particle  size  and/or  shape. 

It  should  be  noted  that  a  wide  range  of  materials  falling  within  the 
class  of  ductile  materials  as  applied  to  dust  erosion  would  be  fortunate 
from  the  viewpoint  of  an  investigator,  since  the  mechanism  of  erosion  of 
brittle  materials  appears  to  be  difficult  to  explain.  This  mechanism  appears 
to  involve  displacement  of  material  by  the  abrasive  particle,  but  in  addition 
the  propagation  and  intersection  of  cracks  makes  an  important  contribution. 

A  detailed  analysis  of  this  latter  mechanism  is  not  available  at  the  present 
time. 
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E.  Concentration  Effect 


Several  tests  were  made  to  determine  the  effect  of  dust  concentra¬ 
tion  upon  the  erosion  rate  of  various  materials.  Intuitively,  it  would  be 
expected  that  the  weight  loss  factor  would  be  independent  of  dust  concen¬ 
tration;  that  is,  an  increase  in  dust  concentration  would  cause  a  propor¬ 
tionate  increase  in  erosion  weight  loss,  and  the  weight  loss  factor  would 
remain  constant. 

Figures  21  and  22  show  the  experimental  results.  It  is  seen  that 
there  is  a  definite  tendency  for  the  weight  loss  factor  to  increase  as  con¬ 
centration  increases.  Although  the  change  in  the  weight  loss  factor  is 
relatively  small,  the  phenomenonis  of  interest,  and  some  possible  explana¬ 
tions  will  be  discussed. 

The  most  obvious  reason  for  the  concentration  effect  observed  is 
that  the  increased  concentration  causes  mutual  interference  between  parti¬ 
cles  such  as  collisions  or  agglomeration  of  several  particles  into  a  large 
low-velocity  mass.  If  this  were  actually  occurring,  however,  the  percent¬ 
age  reduction  in  weight  loss  for  a  given  increase  in  concentration  should  be 
the  same  for  each  material.  From  the  curves,  it  is  seen  that  this  does 
not  hold  true,  when  comparing  the  results  for  aluminum  with  those  for  the 
steel  specimens.  The  percentage  change  in  the  weight  loss  factor  for 
aluminum  is  much  larger  than  for  the  steels. 

As  the  dust  concentration  is  increased  and  all  other  variables  held 
constant,  the  net  effect  upon  the  specimen  is  a  higher  frequency  of  impact 
of  the  dust  particles.  The  cause  for  the  observed  result  therefore  must 
be  related  to  the  time  rate  of  change  of  some  variable  that  affects  erosion 
rate. 


A  possible  explanation  of  the  observed  results  involves  the  work 
hardening  of  the  specimen  material  in  the  neighborhood  of  the  particle 
impact.  As  is  discussed  in  another  section  of  this  report,  the  heat  gener¬ 
ated  by  the  impacting  particle  results  in  rather  high  temperatures  in  the 
specimen  material  near  the  cutting  tip  of  the  particle.  It  seems  reason¬ 
able  to  assume  that  any  work  hardening  of  the  specimen  due  to  the  deforma¬ 
tion  of  the  surface  by  the  impacting  particle  will  decrease  with  time,  due 
to  the  high  temperature  levels.  A  second  particle  impacting  in  the 
deformed  region  immediately  after  the  first  impact  would  encounter  a 
material  of  higher  hardness  than  a  particle  that  impacted  at  some  later 
time.  Thus  a  high  dust  concentration  may  maintain  a  higher  specimen 
surface  hardness  than  a  low  dust  concentration,  with  a  resulting  lower 
rate  of  erosion. 


Weight  Loss  Factor  g:/in 
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A  second  explanation  involves  consideration  of  the  formation  of 
oxides  on  the  specimen  surface.  Cutwater  and  Shaw^,  while  performing 
grinding  tests  in  an  inert  atmosphere,  observed  an  increase  of  grinding 
forces  up  to  25  times  that  experienced  in  an  atmosphere  of  air.  This 
result  was  attributed  to  the  absence  of  an  oxide  film  on  the  freshly  cut  sur¬ 
faces  in  the  inert  atmosphere,  resulting  in  the  rewelding  of  small  chips 
to  the  specimen  surface  as  they  curled  back  onto  the  surface  during  further 
cutting.  This  suggests  an  explanation  for  the  concentration  effect  noted 
during  erosion  tests.  From  the  size  distribution  of  the  0-74  micron 
dust,  it  can  be  calculated  that  there  are  approximately  4.5  X  10^  dust  par¬ 
ticles  over  1,8  microns  in  diameter  per  gram  of  dust.  For  the  test  con¬ 
ditions  of  Figures  21  and  22,  it  may  be  estimated  that  roughly  from  1  to 
100  particles  over  1.8  microns  in  diameter  strike  within  an  area  of  100 
square  microns  each  second,  the  range  of  impact  frequency  being  due  to 
the  variation  in  the  test  concentration.  It  seems  reasonable  to  assume 
that  as  the  impact  frequency  (concentration)  increases,  more  and  more 
often  there  will  be  occasions  in  which  insufficient  time  exists  between 
impacts  for  an  oxide  film  to  form,  resulting  in  minute  particles  of  mate¬ 
rial  rewelding  to  the  surface  with  the  corresponding  decrease  in  weight 
loss. 


The  experimental  data  does  not  permit  a  choice  between  these 
theories..  It  is  interesting  to  note,  however,  that  the  aluminum  alloy, 
which  is  well  known  for  its  ability  to  rapidly  form  an  oxide  film,  exhibits 
a  greater  change  in  weight  loss  factor  per  unit  change  in  concentration  than 
does  steel.  This  could  be  interpreted  to  mean  that  for  a  decrease  in 
impact  frequency,  the  aluminum  is  able  to  form  a  protective  film  more 
rapidly  than  steel,  resulting  in  a  greater  relative  decrease  in  the  alumi¬ 
num  erosion  loss. 

There  is  a  possibility  that  further  study  of  the  observed  concentra¬ 
tion  effect  would  increase  understanding  of  the  mechanism  involved,  and 
the  subsequent  control  of  the  mechanism,  perhaps  by  proper  choice  of 
materials  or  material  coatings,  could  result  in  a  substantial  decrease  in 
erosion  rates. 

F.  Dust  Size  Effect 


Several  series  of  tests  were  conducted  to  define  the  relative  effect 
of  the  different  dust  size  ranges  upon  the  rate  of  erosion,.  A  typical  set  of 
results  is  shown  in  Figure  23,  This  chart  shows  the  variation  in  weight 
loss  factor  for  a  constant  velocity.  From  previous  discussions, 

it  will  be  evident  that  these  results  will  not  accurately  reflect  the  erosion 
capabilities  of  the  dust  at  a  constant  parti  cle  velocity  since  the  larger 
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particles  do  not  attain  airstream  velocity.  However,  with  information 
previously  obtained,  it  is  possible  to  correct  the  test  results  to  show  the 
size  effect  at  a  constant  particle  velocity.  This  is  done  by  using  Equation 
(<21),  rewritten  below 

w  =  + - +  ^  (21) 

Equation  (21)  represents  the  contributions  of  increments  of  size  ranges  of 
dust  to  the  weight  loss  factor.  Five  size  increments  were  used,  the  first 
being  the  size  range  represented  by  the  0-5  micron  dust,  and  the  remaining 
four  were  the  size  range  increments  represented  by  the  difference  between 
the  maximum  particle  diameters  of  the  0-5  and  0-10  micron  dust,  the 
0-10  and  0-43  micron  dust,  and  so  on.  For  each  size  range  of  dust. 
Equation  (21)  was  written.  The  constant  C  was  obtained  from  Figure  3  or 
from  Equation  (5).  The  weight  of  dust  m  in  the  size  increments  was 
obtained  from  Figure  2.  Thus,  for  the  0-5  dust,  Equation  (21)  is  written 
containing  only  one  term 


m^kjC2w2 


For  the  0-10  micron  dust.  Equation  (21)  contains  two  terms,  and  so  on, 

The  weight  loss  factor  w  for  each  case  is  taken  from  the  experimental  data. 
The  five  equations  have  five  unknowns,  the  k  factors.  These  are  obtained 
by  solution  of  the  equations;  the  equations  are  rewritten  letting  the  con¬ 
stants  C] ,  C^,  C3  be  unity,  and  the  corrected  weight  loss  factor  is  obtained. 
This  corrected  weight  loss  factor  represents  the  erosion  for  a  constant 
particle  velocity.  The  corrections  are  shown  in  Figures  23  and  24  as 
dashed  lines.  The  corrections  for  the  0-5  and  0-10  micron  dust  were 
negligible , 

Figure  23  shows  that  the  weight  loss  factor  first  increases  rapidly 
as  the  dust  size  range  increases,  and  then  becomes  fairly  constant  for 
size  ranges  greater  than  0-43  microns.  At  first,  this  seems  to  indicate 
that  the  0-150  micron  dust  is  not  significantly  more  damaging  than  0-43 
micron  dust.  This  statement  is  true,  as  long  as  the  erosion  rate  is  speci¬ 
fied  in  terms  of  unit  weight  of  dust.  That  is,  all  other  variables  held  con¬ 
stant,  a  gram  of  0  -  150  dust  will  cause  approximately  the  same  erosion  as 
a  gram  of  0  -  43  micron  dust.  More  significant,  however,  is  a  plot  of  the 
erosion  rate  on  a  per  impact  basis.  Using  the  data  shown  on  the  dust  size 
distribution  curve,  Figure  2,  and  as suming  particle  sphericity,  the  number 
of  dust  particles  per  gram  of  dust  for  each  size  range  of  dust  may  be 
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FIGURE  23.  WEIGHT  LOSS  FACTOR  VERSUS  DUST 

PARTICLE  SIZE  RANGE,  AIR  VELOCITY 
850  FPS,  DUST  CONCENTRATION  0.  013 
GR/FT3,  peARLITIC  C-  1050  STEEL. 

40"  TEST  IMPACT  ANGLE,  SILICA  FLOUR 


FIGURE  24.  AVERAGE  WEIGHT  LOSS  PER  IMPACT 
VERSUS  DUST  PARTICLE  SIZE  RANGE, 
AIR  VELOCITY  850  FPS,  DUST  CONCEN¬ 
TRATION  0,  0 13  GR/FT^,  PEARLITIC 
C-1050  STEEL  40  “  TEST  IMPACT  ANG.LE, 
SILICA  FLOUR 
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calculated.  Dividing  the  weight  loss  factor  by  the  number  of  particles  per 
gram  of  dust  gives  an  average  weight  loss  per  particle.  Figure  Z4  shows 
a  plot  using  this  erosion  parameter.  The  results  now  show  continuously 
increasing  erosion  with  increasing  size  range,  a  result  intuitively  expected. 

The  theoretical  equations  for  volume  removal  by  a  single  particle 
[Eqs.  (16)  and  (17)]  indicate  that  for  a  particle  with  a  moment  of  inertia 
proportional  to  mr^,  the  volume  removed  is  proportional  to  the  particle 
mass,  or  to  the  diameter  cubed.  It  is  expected  that  tests  using  particles 
of  one  diameter  (instead  of  a  range  of  diameters)  will  show  this  trend,  pro¬ 
viding  the  particles  are  relatively  similar  in  geometry.  It  is  emphasized 
that  the  data  shown  here  are  for  the  specific  size  ranges  of  dust  shown  in 
Figure  2,  and  should  not  be  used  to  predict  erosion  losses  for  a  particular 
dust  particle  diameter,  or  for  dust  distributions  significantly  different  from 
those  used  in  this  study. 

G.  Material  Properties  Effect 

In  Finnie's  equations  for  volume  removal  by  a  single  abrasive  par¬ 
ticle  [  Eqs.  (16)  and  (17)]  .  it  is  s  een  that  the  effect  of  material  properties 
is  included  by  means  of  the  plastic  flow  stress  term,  p.  The  plastic  flow 
stress  of  a  material  is  often  equated  to  the  material  hardness  measured 
by  the  Brinell  method  or  an  equivalent.  Thus,  it  would  be  expected  that 
the  volume  removed  in  an  erosion  test  would  be  inversely  proportional  to 
the  hardness  of  the  materials  tested,  other  variables  held  constant.  Figure 
25  is  a  plot  of  volume  loss  factor  versus  Knoop  hardness  for  several  mate¬ 
rials.  The  general  shape  of  this  curve  was  verified  by  several  tests  at 
different  air  velocities  and  test  impact  angles,  The  data  shown  are  typical 
of  all  these  tests.  The  plot  indicates  that  the  inverse  relationship  could  con¬ 
ceivably  be  valid  for  very  low  material  hardnesses,  but  is  definitely  not  of 
use  throughout  most  of  the  range  of  hardnesses  of  common  engineering 
materials.  (It  should  be  noted  that  this  difficulty  does  not  affect  the  validity 
of  Finnie's  equations,  and  only  shows  that  the  plastic  flow  stress  used  in 
the  equations  is  not  equivalent  to  material  hardness  for  the  materials  tested.) 
Figure  25  indicates  fair  correlation  between  erosion  rate  and  material  hard¬ 
ness,  but  fails  to  explain  the  high  erosion  rate  for  the  martensitic  steel, 
and  the  low  rate  for  the  tempered  steel. 

In  attempting  to  better  correlate  erosion  rates  with  material  pro¬ 
perties,  several  dimensional  analyses  were  conducted.  Both  mechanical 
and  thermal  properties  of  the  materials  were  used  as  variables.  Variables 
used  in  the  analyses  were  material  hardness,  material  density,  particle 
impact  velocity,  material  thermal  conductivity,  material  specific  heat, 


Volume  Loss  Factor 
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FIGURE  2  5.  VOLUME  LOSS  FACTOR  VERSUS  KNOOP  HARDNESS 
NUMBER,  AIRSTREAM  VELOCITY  835  FPS,  0  -  74 
MICRON  SILICA  FLOUR,  DUST  CONCENTRATION 
0.013GR/FT3,  90 “IMPACT  ANGLE 


FIGURE  26.  VOLUME  LOSS  FACTOR  VERSUS  KNOOP  HARDNESS 
NUMBER  X  %  ELONGATION,  AIRSTREAM  VELOCITY 
835  FPS,  0  -  74  MICRON  SILICA  FLOUR,  DUST 
CONCENTRATION  0.  013  GR/FT^,  90’  IMPACT  ANGLE 
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reduction,  of  area  of  material  as  found  in  standard  tensile  test,  modulus  of 
elasticity,  and  others.  Various  dimensional  groups  were  obtained,  but  no 
satisfactory  correlation  was  found. 

From  the  relative  erosion  rates  of  the  various  steel  specimens,  it 
appeared  that  the  strain  energy  required  to  fracture  was  an  important 
variable.  Since  material  hardness  is  a  measure  of  tensile  strength,  the 
product  of  hardness  and  percent  elongation  represents  an  approximation 
of  the  area  under  the  stress-strain  curve,  or  the  strain  energy  to  fracture. 
Figure  26  shows  a  plot  using  this  variable.  The  C-1050  steel  specimens 
now  have  a  reasonable  relation,  but  the  overall  correlation  is  not  greatly 
improved.  Some  reasons  for  the  difficulties  discussed  here  will  now  be 
considered. 

In  the  past,  considerable  work  has  been  done  on  the  mechanisms 
of  material  removal  in  the  grinding  process'^’®,  with  some  very  interesting 
results.  Comparison  of  the  grinding  process  with  the  erosion  of  materials 
by  small  particles  shows  considerable  similarity  between  the  two.  The 
depth  of  cut  of  the  abrasive  grains  in  the  material  is  of  the  same  order  of 
magnitude  and  the  speed  of  cutting  is  also  similar.  For  the  grinding  pro¬ 
cess,  Backer,  Marshall,  and  Shaw®  have  demonstrated  that  the  specific 
energy  of  material  removal  (the  energy  required  to  remove  a  unit  volume 
of  material)  is  extremely  high,  and  corresponds  to  a  material  strength 
equal  to  the  theoretical  material  strength.  This  theoretical  strength  may 
be  calculated  by  considering  the  attractive  and  repulsive  forces  of  atoms 
in  adjacent  rows  of  a  perfect  crystal  lattice.,  The  theoretical  strength  of 
a  given  material  can  be  shown  to  be  approximately  equal  to  G/ Zir,  where 
G  is  the  shear  modulus  of  the  material.  The  theoretical  strength  is  there¬ 
fore  not  affected  by  changes  in  hardness  or  other  properties  that  do  not 
involve  a  change  in  G,  For  steel,  the  theoretical  strength  is  about  1  .8 
X  10^  psi,  compared  to  values  of  about  50, 000  psi  found  in  tensile  tests 
for  mild  ste<'l.  It  is  thought  that  the  difference  in  these  two  values  may 
be  due  to  the  imperfections  (grain  boundaries,  crystal  defects,  and  impuri¬ 
ties)  actually  present  in  the  material.  However,  in  the  grinding  process, 
the  depth  of  cut  is  of  the  same  order  of  magnitude  as  the  mean  distance 
between  imperfections,  and  the  material  is  able  to  demonstrate  its  theoret¬ 
ical  strength.  This  effect  is  shewn  in  Figure  27  for  the  grinding  process. 
These  results  indicate  that  for  depths  of  cut  less  than  about  28  microinches, 
the  specific  energy  remains  constant  at  its  theoretical  value. 

For  sev-rai  dust  erosion  tests,  the  depth  of  erosion  pits  on  the 
specimen  surf;  -es  were  measured  with  a  Talysurf  surface  measuring 
instrument,  i  he  aver^.ge  pit  depths  for  steel  specimens  ranged  from  about 
20  microinches  for  the  Icwer  air  velccities  to  50  microinches  for  air 


Depth  of  Cut  -  Microinches 

FIGURE  27.  SPECIFIC  ENERGY  OF  VOLUME  REMOVAL 
VERSUS  DEPTH  OF  CUT  FOR  GRINDING 
TESTS  *  BACKER,  MARSHALL  AND  SHAW. 
REFERENCE  8 
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velocities  of  850  fps.  It  seems  likely  that  a  significant  percentage  of  the 
total  number  of  impacts  result  in  pit  depths  much  small  than  this,  since 
the  smaller  pits  are  obscured  by  the  larger  in  the  surface  analysis.  It 
would  therefore  be  expected  that  the  specific  energies  involved  would  be 
high,  approaching  the  theoretical  limit.  This  is  indicated  in  the  present 
study.  The  specific  energy  for  several  tests  was  calculated  by  computing 
the  particle  kinetic  energy  used  in  removing  unit  volume  of  material.  For 
steel,  the  values  obtained  ranged  from  about  3  X  10^  in. -lb  per  cu  in.  to 
9  X  10^  in. -lb  per  cu  in.  The  fact  that  these  values  are  higher  than  those 
obtained  in  grinding  tests  is  reasonable  when  it  is  considered  that  many  of 
the  paiLicles  are  ineffective  ii'i  remov  uig  nr  ate  rial  due  to  unfavorable 
impact  angles,  and  that  many  of  the  particles  rebound  from  the  specimen 
retaining  a  portion  of  their  original  kinetic  energy. 

The  effect  of  these  considerations  upon  the  erosion  mechanism 
may  be  visualized  as  follows:  The  smaller  dust  particles,  when  able  to 
remove  material  from  the  specimen,  are  forced  to  do  so  at  stresses  near 
the  theoretical  strength  of  the  material.  The  larger  particles  are  able  to 
cut  the  material  at  lower  stresses,  dependent  on  the  depth  of  cut  and  the 
properties  of  the  material.  To  obtain  a  better  correlation  between  mate¬ 
rial  properties  and  erosion  rate,  it  would  seem  that  the  most  straightfor¬ 
ward  approach  would  be  a  series  of  tests  using  particles  of  essentially 
the  same  size  for  each  tfjst,  rather  than  a  wide  range  of  sizes  as  was  done 
in  this  investigation.  This  would  enable  the  investigator  to  separate  the 
erosion  occur  ring  at  t!ie  theoretical  limit  due  to  small  particles  from  that 
occurring  at  lesser  stresses  due  to  the  larger  particles,  A  more  favorable 
condition  for  determining  the  material  properties  controlling  the  latter  type 
of  erosion  would  then  exist. 

The  results  of  material  property  studies  as  applied  to  d\ist  erosion 
will  be  influenced  by  the  temperatures  generated  by  the  cutting  particle. 
Outwater  and  Shaw"^  have  shown  that  for  grinding  steel,  temperatures  over 
^000°  F  are  probable  in  the  neighborhood  of  the  cutting  edge  of  the  particle. 
During  dust  erosion  tests  on  a  martensitic  C-]050  steel  specimen,  sparks 
were  observed  downstream  from  the  specimen.  These  sparks  seemed  to 
originate  about  six  inches  from  the  specimen  and  caused  a  brief  streak  of 
light.  The  presence  of  these  sparks  indicates  a  high  material  temperature 
due  to  the  friction  forces  of  cutting.  It  would  be  expected  therefore,  that 
the  high  temperature  properties  of  the  material  will  be  significant  in 
determining  the  characteristics  of  materials  with  respect  to  erosion  by 
dust. 


Outwater  and  Shaw^  indicate  that  the  mean  temperature  attained 
during  the  grinding  process  is  a  function  of  the  shear  energy  of  the  cutting 
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process,  the  specific  heat  of  the  material,  the  geometry  of  the  abrasive 
particle- specimen  system,  and  the  depth  of  cut.  Applied  to  the  erosion 
process,  it  may  be  expected  the  temperature  will  not  only  vary  from  mate¬ 
rial  to  material,  but  also  will  vary  for  a  given  material  from  impact  to 
impact  due  to  different  depths  of  cut  of  the  various  sizes  of  particles. 

To  summarize  the  results  discussed  in  this  section  of  the  report: 

(1)  For  the  materials  tested,  material  hardness  does  not  ade¬ 
quately  define  the  plastic  flow  stress  used  in  Finnie's  equation, 

(2)  For  erosion  of  materials  by  air-borne  dust,  the  material 
theoretical  strength  is  believed  to  have  a  considerable  effect 
on  the  erosion  rate, 

(3)  A  fair  degree  of  correlation  exists  between  erosion  rate  and 
material  fracture  energy. 

(4)  Better  correlation  of  material  properties  versus  erosion  rates 
will  probably  require  knowledge  of  the  high  temperature  pro¬ 
perties  of  the  material. 
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IV.  CONCLUSIONS  AND  RECOMMliJNDATIONS 


From  the  results  cf  the  study  presented  in  this  report,  it  is  concluded 

that: 

(1)  The  erosion  by  air-borne  dust  of  the  materials  considered  is 

a  process  of  plastic  displacement  of  material  by  the  impacting 
dust  particle,  and  the  effects  of  dust  particle  velocity  and  angle 
of  impact  upon  erosion  rates  may  be  adequately  described  by 
equations  defining  the  trajectory  of  the  particle  through  the 
material . 

(2)  Contrary  to  previous  acceptance  the  concentration  of  dust  in 
the  airstream  has  a  significant  effect  upon  the  weight  loss 
factor. 

(3)  Erosion  rates  are  a  strong  function  of  dust  particle  size,  the 
correlation  being  best  described  by  consideration  of  the  average 
erosion  loss  per  particle  impact. 

(4)  The  erosion  rate  of  materials  is  probably  affected  by  the 
theoretical  material  strength  for  the  smaller  range  of  particle 
sizes.  A  detailed  examination  of  the  pertinent  material  pro¬ 
perties  applied  to  dust  erosion  should  consider  this,  and  also 
the  high  temperature  properties  of  the  material. 

It  is  recommended  that  further  work  be  done  to  increase  understanding 
of  the  relationship  betwet  n  erosion  rate,  particle  size,  and  material  properties. 
It  is  believed  that  this  can  be  carried  out  most  efficiently  by  experimental 
work  patterned  after  the  present  study,  but  with  the  addition  of  equipment  to 
separate  dust  into  small  ranges  of  size.  This  control  ever  dust  size  would 
eliminate  many  of  the  difficulties  encountered  in  the  present  investigation, 
and  would  allow  the  investigator  to  control  more  closely  impact  velocity, 
impact  angle,  and  particle,  depth  of  cut.  Such  an  investigation,  it  is  believed, 
would  provide  significant  information  on  the  effect  of  material  properties 
upon  erosion  rate. 

For  application  to  the  dust  erosion  cf  gas  turbines,  it  is  considered 
essential  that  further  experimentation  be  done  at  elevated  specimen  tempera¬ 
tures,  and  the  effect  of  specimen  temperature  upon  erosion  rate  be  defined. 

In  addition,  the  relative  abr asi^'enes  s  of  different  types  of  natural  dust  needs 


to  be  determined  for  purposes  of  practical  design.  From  the  results  of 
this  investigation,  it  is  seen  that  to  completely  define  the  properties  of  a 
dust  as  pertains  to  erosion,  three  relations  must  be  given.  First,  the 
erosion  rate  as  a  function  of  particle  angle  of  impact  for  the  dust  will  pro¬ 
vide  information  as  to  the  magnitude  of  the  angle  for  maximum  erosion 
rate.  Second,  the  relationship  between  erosion  rate  and  particle  size 
must  be  defined.  This  will  require  separation  of  the  dust  sample  into 
narrow  size  ranges.  Third,  from  the  results  of  the  first  two  tests,  the 
magnitude  of  the  weight  loss  factor,  or  a  similar  parameter,  must  be  deter 
mined  relative  to  a  standard  dust  eroding  a  standard  specimen.  It  is  pos¬ 
sible  that  many  natural  dusts  will  have  common  values  for  one  or  more  of 
these  relationships;  however, this  should  be  verified  experimentally. 
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APPENDIX  A 

DERIVATION  OF  PARTICLE  DISTANCE- 
VELOCITY  RELATIONSHIP 


For  a  particle  being  accelerated  downward  by  a  fluid  stream: 

^  .  f  P  -  Po\  dV 

F+mg(— )  -m— 

where  the  drag  force  F  is 

F  =  ^  Po  (W  -  V)^  ACr 

The  particle  Reynold's  number  is  given  by 


(Al) 


(A2) 


Po'i 

R  =  —  (w  -  V) 

The  initial  particle  Reynold's  number  is 


Define  the  term  p  as 


Zmgpod  ^  p  -  Po  \ 


and  let  the  drag  coefficient  be 


(A3) 


(A4) 


(A5) 


>  40 

Cr  -  0.4  +  — 


(A6) 


Combining  Equations  (a1)  through  (A6)  results  in 

2md  -dR 
dt  = — = - - - 


(A7) 
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Integrating  (A7)  from  Rq  to  R  letting 

S  =  'n/'402TT7^  (A8) 

results  in 


t  = 


2md 


log 


.  8R  +  40  -  S  .SRq  +  40  +  S  \ 
.  8R  +  40  +  S  8Ro  +  40  -  S  / 


(A9) 


For  the  range  of  air  properties  and  particle  sizes  used  in  this  report, 
Equation  (A8)  reduces  to 


S  »  40 


(AlO) 


Assuming  that  the  particles  are  spheres 


m 


p  and 


A 


(All) 


Equation  (A9)  may  be  rearranged  to  give  V  explicitly,  and  integrated  to 
give  the  distance  x  traveled  by  the  particle  as  a  function  of  time  t.  Sub¬ 
stituting  Equations  (A9),  (AlO),  and  (All)  into  the  equation  for  x,  and 
rearranging 


Pq  X 
P  d 


^o  r  R  +  100 
30  Rq  +  100 


+  3  j  log 


R  +  100 
Ro  +  100 


(A12) 
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APPENDIX  B 

ESTIMATION  OF  TRUE  IMPACT  ANGLE 


A  sketch  of  the  specimen  in  the  airstream  is  shown  below, 


The  airstream  is  shown  schematically,  and  a  typical  particle  trajectory 
is  represented.  The  number  of  particles  impacting  the  specimen  at  any 
location  L  per  unit  specimen  area  per  unit  time  is 

N  =  CV  sin  a  (Bl) 

where  C  and  V  are  the  local  particle  concentration  and  velocity  respectively. 
C  and  V  depend  on  the  conditions  of  air  flow  at  tne  specimen,  but  for  a  small 
specimen  with  a  corresponding  small  airstream  disturbance,  C  and  V  may 
be  approximated  by  a  constant.  Then 

sin  a  =  kN  (Bil) 

where  k  is  a  constant,  From  consideration  of  the  air  flow  and  the  resulting 
particle  trajectories,  it  is  evident  that  at  some  location  on  the  specimen, 
there  will  exist  a  stagnation  point  where  the  impact  angle  a  is  equal  to  the 
specimen  angle  0.  Let  this  location  be  designated  by  Lq,  At  Lq,  then,  N 
will  be  a  maximum,  since  there  is  no  dilution  of  the  particle  concentration 
by  the  bending  of  the  airstream.  This  may  be  written 

when  L  =  Lq.  a  -  0,  N  =  Nmax  (B3) 

Substituting  (B3)  into  (B2)  results  in 
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,  sin  0 
k  =  — - 

A  V 


Su\stituting  for  k  in  (Bi)  gives 


sin  a  = 


sin  9 


